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Abstract 

The circulation around any closed loop is a Lagrangian invariant for classical, 
smooth solutions of the incompressible Euler equations in any number of space 
dimensions. However, singular solutions relevant to turbulent flows need not pre- 
serve the classical integrals of motion. Here we generalize the Kelvin theorem on 
conservation of circulations to distributional solutions of Euler and give necessary 
conditions for the anomalous dissipation of circulations. We discuss the important 
role of Kelvin's theorem in turbulent vortex-stretching dynamics and conjecture a 
version of the theorem which may apply to suitable singular solutions. 
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1 Introduction 



In a monumental paper, iHelmholtzl (|l858h formulated the fundamental laws 



of vortex motion for incompressible fluids. These include the statements, in 
three space dimensions, that vortex lines are material lin es and that th e flux 
within any vortex tube is a Lagrangian invariant. Lord Kelvin ( 1869^ gave 



an elegant alternative formulation of these laws in terms of the conservation 
of circulation, for any closed loop advected by an ideal fluid. This theorem 
is equally valid in any space dimension. However, it is o nly rigorously pr oved 
for sufficiently smooth solutions. As was pointed out bv lOnsager (Ell), the 



classical conservation laws need not be valid for singular solutions of Euler 
equations. In particular, breakdown of the energy conservation law can ac- 
count for the anomalous dissipation of energy observed in turbulent fluids at 
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high Reynolds numbers. From a physical point of view, this breakdown of 
energy conservation corresponds to the turbu lent energy cascade and a flux 



of ene r gy to arbitrarily small sc ales. See also Evink f 1994V Constantin et al. 
f)l994l ). lDuchon fc Robertl (fcOOOh . 



These considerations make it very natural to inquire whether Kelvin's theorem 
will remain valid for singular solutions of the Euler equations. This question 
asumes some importance since the conservation of circulations was argued by 
Tavlor ( 1938f ) to play a key role in the enhanced production of dissipation in 
turbulent fluids, by the process of vortex line-stretching. Despite its plausibil- 
ity, the validity of Taylor's argument is far from clear. It is not obvious a priori 
why there should not be anomalous dissipation of the circulation invariants, 
corresponding to a turbulent "flux of circulations" from large to small scales. 

In this paper, we examine these questions and establish a few relevant rig- 
orous results. In section 2 we briefly review the classical Kelvin-Helmholtz 
theorem and its role in turbulence dynamics. In section 3 we prove a theorem 
on conservation of circulations for singu l ar sol utions of incompressible Euler 
equations, analogous to that of lOnsager (1949) for conservation of energy. In 
section 4 we discuss difficulties in formulating the Kelvin theorem for singular 
solutions, due to the breakdown in uniqueness of Lagrangian trajectories, and 
conjecture a statistical version of circulation-conservation which may apply. 



2 The Classical Kelvin Theorem 



The velocity field u(x, t) solving the incompressible Navier-Stokes equation 
a t u + (u'V)u=-Vp+i/Au, V-u = (1) 



with x G A C R d , for any integer d > 2, satisfies the Kelvin-Helmholtz 
theorem in the following sense: For any closed, rectifiable loop C C A at an 
initial time to, the circulation T(C,t) = § C t t \ u(t)'dx satisfies 

^T{C,t) = u f Au(f).dx, (2) 

C(t) 



sec 



where C(t) is the loop advected by the fluid velocity, at time t. E.g., 
Saffman ( 19921 ). section §1.6, for the standard derivation. It is worth observing 



that the Kelvin theorem for all loops C is formally equivalent to the Navier- 
Stokes equation (1). Indeed, if u(x, t) is a smooth spacetime velocity field, 
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diverence-free at all times t, then equation (2) implies that 

j> [D t u(t) - v A u(t)] -dx = (3) 

c 

for all loops C at every time t. Here D t u = d t u + (u- V)u is the Lagrangian 
time-derivative and the equation (3) is derived by applying (2) to the pre- 
image of the loop C at initial time to. By Stokes theorem, equation (3) can 
hold for all loops C C A if and only if there exists a pressure-field p(pc, t) 
such that the Navier-Stokes equation (1) holds locally and also globally, if the 
domain A is simply connected. 

In the inviscid limit v — > 0, the circulation is formally conserved for any initial 
loop C. The fluid equations in this limit, the incompressible Euler equations, 
are the equations of motion of a classical Hamiltonian system. They can be 
derived by the Hamilton-Maupertuis principle from the action functional 

1 tf 

S[x] = -JdtJdeL |x(a,t)| 2 (4) 

to A 



with the pressure field p(x, t) a Lagrange multiplier to enforce the incom- 
pressibility constraint. Here x(a, t) is the Lagrangian flow ma p which satisfies 
x(a, t) = u(x(a, t),t) with initial condition x(a, t ) = a. See Salmon ( 19881) 
for a review. This variational principle yields the fluid equations in a La- 
grangian formulation, as x(a, t) = — Vp(x(a, t), t). The Eulerian formulation 
(1) (with v = 0) is obtained by performing variations in the inverse map 
a(x, t), or "back-to-labels map", with fixed particle positions x. This Hamil- 
tonian system has an infinite-dimensional gauge symmetry group consisting of 
all volume-preserving diffeomorphisms of A, which corresponds to all smooth 
choices of initial fluid particle labels. In this framework, the conservation of 
the circulations for all closed loops C emerges as a consequence of Noether's 
theorem for the particle- relabelling symmetry. See Salmon (1988), Section 4. 



The Kelvin theorem has many profound consequences for the dynamics of 
incompressible fluids. We just note here a w ell-known deduction for three- 



I. 



dimensional turbulence by G 
vortex lines — or, for that matter, any material lines 



Tavlor (1938). It is reasonable to assume that 
will tend to elongate un- 
der chaotic advection by a turbulent velocity field. Incompressibility requires 
that the cross-sectional area of a vortex tube formed by such lines will shrink 
with time. But, in that case, the Kelvin- He lmholtz theor em implies that the 
vorticity magnitude of the tube must grow. Tavlor ( 19381 ) observed that this 
process of vortex line-stretching provides an intrinsic mechanism for ampli- 
fication of the net viscous dissipation v f A dx.\u>\ 2 , where u) = Vxu. More 
recently, regularizations of the Navier-Stokes equation have been proposed as 
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model equations for large-scale turbulence by iHolm et al.l (|1998T ) ; iFoias et al. 
f)2001l l. motivated by requiring that a Kelvin circulation theorem be preserved. 



3 A Generalized Theorem on Conservation of Circulation 



This line of reasoning of Tavlor ( 1938j ) is quite delicate. It assumes certain 



properties (material advection of vortex lines and conservation of circulations) 
that can hold strictly only in the v — > limit. However, conclusions are drawn 
about the limiting behavior of the energy dissipation rat e, which directl y in- 
volves the kinematic viscosity v\ Furthermore, as noted bv lOnsager (|l949h . the 



solutions of the Navier-Stokes equation are not expected to remain smooth in 
the inviscid limit. Thus, the righthand side of equation (2) does not necessarily 
vanish as v — > 0. Physically, there may be a dissipative anomaly for the con- 
servation of circulations. If so, then the validity of Taylor's vortex-stretching 
mechanism for turbulent energy dissipation is open to serious question. 

The formulation of turbulent conservation of circulation by a zero-viscosity 
limit is physically natural, but not the most convenient either for numerical 
tests or for rigorous mathematical analysis 1 . We shall instead consider directly 
the singular (or distributional) solutions u G L 2 ([0, T], A) of the incompressible 
Euler equations, with v = 0. Let = Gg * u denote the low-pass filtered 
velocity at length-scale £, where Gi(r) = £~ d G(r/£) is a smooth filter kernel. 
Then satisfies the following equation (in the sense of distributions in time): 

d t U£ + (u r V)u £ = - Vp e + f e , (5) 



where p e is the filtered pressure and where = — V-r^ is the subgrid force, 
i.e. minus the divergence of the stress-tensor Tg = (uu)^ — u^u^. Let us choose 
a rectifiable closed loop C in space. We define Ce(t) as the loop C advected by 
the filtered velocity u^. This definition makes sense, since the filtered velocity 
Uf is Lipschitz in s pace, and corresponding flow maps exist and are unique 
(DiP erna fe Lions! (|l989i )). We define a "large-scale circulation" with initial 
loop C as the line-integral Ti(C,t) = ■ dx. The same calculation 

that establishes the Kelvin circulation theorem for smooth solutions of Euler 



1 For example, the Kelvin theorem in the form of equation (2) is not proved to be 
valid f or the global solutions of the Navier-Stokes equation (1) constructed bv lLerav 
( 19341 ). The difficulty here is that the Leray regularization of (1) does not preserve 
the Kelvin theore m, while alternative regularizations wh ich do, such as t hat of 
Foias et all (|200lh . in turn modify the energy balance. See IConstantirJ fcooah . 
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equations gives that 



r e (C, t) - T e (C, t ) = J dr j f,(r) ■ dx. 

to C e (r) 



(6) 



Thus, the line-integral of i? on the RHS represents a "flux" to subgrid modes 
at length-scales < i of circulation on the loop Cgir). This motivates the defi- 
nition, for any loop C and filter length £, 



K e (C,t) = - j f e (t)-dx 



(7) 



C(t) 



so that (in generalized sense) (d/dt)Ti(C,t) = —K^(C,t). 
We now prove the following: 

Theorem: Let( p be the pth-order scaling exponent of the velocity, in the sense 
that it is the maximal value such that 



i-|^x|5u(r;x)r = 0(|rH, 



for all |r| < r , where <5u(r; x) = u(x + r) — u(x). Then for any smooth loop 
C C A 

K t (C) = -ff t -dx 
c 

satisfies lim^ Kg(C) = if ( p > (d — 1) + (p/2) for any p > 2. 

The special case of this result for p = oo states that the "circulation flux" will 
go to zero as £ — > if the smallest v elocity Holder e xponent h min is > 1/2. This 
is an exact analogue of the result of Onsager ( 1249) for vanishing of energy flux 
when h m i n > 1/3. One can see that it is even easier for circulation-conservation 
to be anomalous than for energy-conservation. 



Proof: Our argument is close to that given by IConstantin et al.l ([19941) for the 
Onsager theorem. The following identity for the subgrid force is easily verified: 



/.(x)=/A( aj «)(r)Mr;x)^(r;x) 

-/A M (r| k (r;x)/A' C M h ,(r';x) 



We omit here all subscripts i for convenience. By this identity, 
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tf'dx= J d d Y(djG){r) 1 6ui(T)6uj(T)da 
c ' Ic 

- J d d r{djG){r) Jd d r'G{r') j> ' 5ui{v)5uj{T')d. 



Thus, 



f • r/.x < J d d r\VG(r)\ j |5u(r)| 2 rfs 

+ J d d r\VG(r)\ J dVG(r') j> \5vl{t)\ \5u{r')\ds 



where s denotes arclength along the curve C. By normalization J d d r' G(r') 
and the inequality |5u(r)| |<5u(r')| < | [|<5u(r)| 2 + |5u(r')| 2 ] , this becomes 



f • dx 



<^|^r|VG(r)| f\5u(r)\ 2 ds 
ic 

+ ^ y d d r|VG(r)| y A'G(r') ^ |<5u(r')| 2 ds 



We now use the Holder inequality to derive the bound 
j> \5u(r)\ 2 ds < [L(C)] (p - 2)/p j \5u(r)\ p ds 



(9) 



c 



\c 



for any p > 2, where L(C) is the length of the curve C. The condition on 
the scaling exponent in the statement of the theorem can rephrased as the 
condition that u belong to the Besov space (A) with o v = ( p /p. (More 
properly, we should replace a v by <j v — e for any small e > 0.) Standard 
trace theorems then imply that the restriction of u to the subman i fold C 
of codimension d-1 must satisfy u| G G See iTriebell ()l983lh 

Theorem 2.7.2. Together with inequality (9), this implies that § c |<5u(r)| 2 <is = 
O ( I x-l 2 tCp i)]/jp ] 2 _ jf b ounc [ j s substituted into estimate (8) for circulation- 



2 Only the case p = oo rigorously follows from standard trace theorems. The prob- 
lem is that the intrinsic Besov space norms on the submanifold C measure only 
the increments between points both on C. However, existing trace theorems im- 
ply that every element / G Bp^oo(C), a' p = a p — (d — l)/p is the restriction to 
C of some element / G Bp%{K). The result we need follows if the semi-norm 
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flux, it yields 



f • G?X 



C 



q ^2[C P -(d-l)]/p-l^ _ 



We see that the latter goes to zero as t — > 0, if ( p > p/2 + (d— 1). □ 



As an application, consider a velocity field that is Lipschitz regular, so that 
( p = p for all p > 1. In that case, it suffices to take p > 2(d — 1) in order to 
show that the circulation is conserved (in the sense that the flux vanishes for 
t — > 0.) This result applies to the 2D enstrophy cascade, since it is ex pected 
there that ^ p = p, with only logarithmic corrections, for all p > 2. See lEvink 
(2000). Thus, we expect that the Kelvin theorem holds in a strong sense - 
for individual realizations — in the 2D enstrophy cascade. However, in the 3D 
energy cascade the conditions of the theorem are not expected to be satisfied. 



4 The Role of Kelvin's Theorem for Singular Solutions 



Even assuming that the assumptions of our theorem are met, there are addi- 
tional difficulties in justifying constancy of the circulation invariants. Vanish- 
ing of the circulation flux for loops of finite length, as established in our theo- 
rem, is not sufficient. In the first place, the material loop C^{t) is not expected 
to remain rectifiable as i — > 0, but instead to becom e a fractal curve C (t) with 
Hausd orff dimension > 1 for any positive time t ( Sreenivasan fc Meneveail 
( 19861 )). Thus, we cannot immediately infer that the RHS of equation (6) 
vanishes as i — > 0, nor even make sense of the contour integral in that 
limit. A possible approach here is to transform the RHS to label-space, as 
§ c if (x>(t), r) • dx.e(r) where the map satisfies x>(a, r) = u^(xg(a, r), r). This 
can make sense as a Stiel ties integral o n the lo o p C i n label-space for Holder 



tie lo o p ( 

continuous maps (e.g. see lYound (|l936t )) IZahlel (fl998f )). 



However, there is a much more serious problem in formulating Kelvin's the- 
orem for singular Euler solutions: it is not clear that material loops exist! 
Recent work on an idealized turbulence problem — the Kraichnan model of 
random advection — has shown that Lagrangian particle trajectories x(t), x'(t) 
can explosively separate even when x = Xq initially, if the advecting v elocit y 
field is only Holder continuous and not Lipschitz. See Bernard et al.l (^998). 
Mathematically, this is a consequence of the non-uniqueness of solutions to the 
initial- value proble m, while, physicall y, it corresponds to the two-partic le tur- 
bulent diffusion of iRIchardsonl (|l92fih . iLe Jan fc Raimondl ([20021 12004 ) have 



SUPj r |< p 



\r\> 



I C ds\Sf(rW 



i/p 



is equivalent to the standard Besov semi- 



norm on B P %{C). Here r ranges over a ball of radius p inside A 
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rigorously proved that there is a random process of Lagrangian particle paths 
x(t) in the Kraichnan model for a fixed realization of the advecting veloc- 
ity and a fixed initial par ticle position. This phenomenon has been termed 



spontaneous stochasticity (IChaves et al. 



alized flow" was proposed bv iBrenierl ( 1989^ for the problem of minimizing 



A similar notion of "gener- 



the action (4). In his formulation, the action is generalized to a functional 
S[P] = \f P(dx) ft* dt |x(t)| 2 , where P is a probability measure on path- 
space, and he showed that minimizers always exist in this framework. Unfor- 
tunately, this notion does not permit one to define the concept of material 
lines and surfaces for ideal flow. A more natural generalization of the classical 
action would be of the form 

*/ 

S[P] = -J P{dx) J dt J da|x(a,t)| 2 (10) 

to A 



where P is now a probability measure on time-histories of measure-preserving 
maps 3 . For any realization of such a random process and for any initial curve 
C the advected object C(t) = x(C, t) is well-defined and remains a (random) 
curve for all time t, if the maps are continuous in space. 

Let us assume for the moment that the (very nontrivial) problem can be solved 
to construct such a generalized flow x(a, t), or stochastic process in the space 
of volume-preserving maps, which is hopefully a.s. Holder continuous in space 
so that material loops C(t) exist as random, fractal curves. We would like to 



3 The Kraichnan model might also benefit from a formulation in terms of maps. 
Formally, a group of Markov transition operators Sf t , can be defined on spaces 
of functionals of maps, with a fixed realization of the velocity u, via a Krylov- 
Veretennikov expansion: 

Sl V = I](-l) n J dt! J dSi! Jdt 2 Jda 2 ---J dtnJdSin 

n=0 t' A f A t' A 

u(x(ai),ii).— 6 —) eft-*^ (u(x(a 2 ),i 2 ).— S j— ) e^Ko 



ox(a n )/ 



Cf. ILe Jan fc Raimondl hm± hooi ). Here the time-integrals should be denned in 



the Ito sense with respect to the white-noise velocity field u(x, t) and Co is formally 
the infinitesimal generator of a diffusion process on the space of maps, given by 
A) = § Ja da Ja da'Aj(x(a) - x(a')) 

5xi(J)6xj(a') ' Gaussian random velocity has 
covariance (ui(x,t)uj(x' ,tf)) = Dij(x — x')5(t — t'). It would be very interesting to 
give rigorous meaning to this expansion, especially for the case where the advecting 
velocity is only Holder continuous but not Lipschitz in space 
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present some plausibility arguments in favor of the conjecture that circulations 
shall be conserved in a statistical sense. More precisely, we expect that the 
circulations T(C,t) for any initial smooth loop C shall be martingales of the 
generalized flow: 

E[T(C,t)\T(C,r),T<f\ =T{C,t'), for t > t'. (11) 



Here E[-} denotes the expectation over the ensemble of random Lagrangian 
paths and we have conditioned on the past circulation history {T(C, r), r < t'}. 
Heuristically, 

(d/dt)E [r(C, t) |r(C, r), t < t'\ = - lim E [K e {C, t) \T{C, r), r < t'} . (12) 



Note that the conditioning event involves scales of the order of the radius of 
gyration of the loops C(r), t < t', while the circulation- flux involves velocity- 
increments over separation lengths I — *■ 0. Therefore, we expect that Kol- 
mogorov's idea of small-scale homogeneity (and isotropy) will apply Note, 
however, that the homogeneous average of the subgrid force ie is zero, because 
it is the divergence of the stress tensor. From another point of view, the sub- 
grid force will become increasingly irregular for i <C R(t) (the radius of the 
loop C(t)) and the sign of the integrand f^(3Q(s, £))oc^(s, t) will oscillate more 
rapidly as a function of the arclength s. Thus, cancellations will occur. For 
these reasons, we expect that the limit on the RHS of (12) shall vanish, im- 
plying (11). Another formal argument can be given by applying the Noether 
theorem to the generalized action (10) and using the fact that a global mini- 
mizer must also minimize the action for the time segment [t', tf]. On the other 
hand, based upon our earlier theorem, it is very unlikely that circulation-flux 
will vanish as I — > in every realization, without any averaging. 

In this section we have clearly indulged in some speculative thinking, but we 
hopefully have also succeeded in outlining the various difficulties in properly 
formulating Kelvin's the orem for turb ulent solutions of the Euler equations. 
Our own view is that the iTavlo r (1938) mechanism of vortex line-stretching is 
the underlying cause of enhanced dissipation in three-dimensional turbulence 
asymptotically at high Reynolds numbers. However, much work remains to 
elucidate the details of the subtle dynamics involved. 
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